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r**- \ A local void in the globally Friedmann-Robertson- Walker (FRW) 

cosmological model is studied. The inhomogeneity is described using the 
Lemaitre-Tolman-Bondi (LTB) solution with the spherically symmetric matter 
distribution based on the faint galaxies number counts. We investigate the effects 
this has on the measurement of the Hubble constant and the redshift-luminosity 
distance relation for moderately and very distant objects (z ~ 0.1 and more). 
The results, while fully compatible with cosmological observations, indicate that 
if we happened to live in such a void, but insisted on interpreting cosmological 



observations through the FRW model, we could get a few unexpected results. 
For example the Hubble constant measuremei 
the separation of the source and the observer. 
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For example the Hubble constant measurement could give results depending on 
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1. INTRODUCTION 

It seems, particularly after the introduction of the inflationary paradigm (Guth 1981), 
that the isotropic and homogeneous Friedmann-Robertson- Walker (FRW) cosmological 
models are best suited for the description of the global structure and the evolution of 
the universe. However, a similar statement is not necessarily true when cosmologically 
moderate scales are thought of. 

There exists direct observational evidence in favour of the large scale isotropy of 
the observed universe, namely, the COBE data confirming a high degree of isotropy of 
the cosmic microwave background radiation (CMBR) (Mather et al. 1990; Hogan 1990; 
Smoot et al. 1993; Wright et al. 1993). However, the observational basis for the other 
standard assumption made in the FRW cosmology, the homogeneity, is weaker. There 
exists observational evidence in favour of larger and larger structures (e.g. Lynden-Bell et 
al. 1988; Clowes and Campusano 1991; Geller and Huchra 1989). 

We think that there exists sufficient observational evidence (discussed later in this 
paper) to support a conjecture that we may live in a relatively large underdense region 
embedded in a globally FRW universe. Exploring observational properties of such a model 
is the aim of the present paper. 

The simple model presented in this paper is restricted by the two physical demands 
we impose on it. Firstly, to be consistent with the CMBR isotropy it has to be spherically 
symmetric. Secondly, the model should be very similar to an FRW one at the beginning of 
the expansion, but become observably different at later times. In this manner, we could 
retain the accomplishments of the FRW cosmology in dealing with early epochs, while 
gaining new freedom in modelling the more recent universe. 

The work on modelling voids of the Lemaitre-Tolman-Bondi type □ in the expanding 
FRW universe has been extensive. Main results and references will be briefly discussed later. 
(An excellent review of the inhomogeneous cosmology exists: Krasihski 1994.) Nevertheless, 
the observable consequences of such a model have seldom been studied (e.g. Paczyhski and 
Piran 1990; Moffat and Tatarski 1992). 

In the following section, we briefly describe the LTB model. Section |3] consists of a 
concise review of the main results in modelling LTB voids, a discussion of the observational 
background and the simple model of a local void presented here. The description of our 



■""A cosmological solution spherically symmetric about one point was first proposed by Lemaitre (Lemaitre 
1931; 1933). However, it is usually called the Tolman-Bondi solution (Tolman 1934, Bondi 1947). 
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results of numerical calculations is contained in Section [I] The closing section contains a 
discussion and conclusions. 

Throughout this paper we use units in which G — c — 1, unless stated otherwise. 
Moreover, we choose the cosmo logical constant A = 0. 



2. MODEL 

First, for the sake of notational clarity, let us recall the FRW line element: 

dr 2 



ds 2 = dt 2 -a 2 (t) 



1 — kr* 



+ r W 



(1) 



with dQ 2 = d6 2 + sin' 



Now, let us consider a Lemaitre-Tolman-Bondi (Tolman 1934, Bondi 1947) model 
for a spherically symmetric inhomogeneous universe filled with dust. The line element in 
comoving coordinates can be written as: 

ds 2 = dt 2 - R' 2 (t, r)f- 2 dr 2 - R 2 (t, r)dVt 2 , (2) 

where / is an arbitrary function of r only, and the field equations demand that R(t, r) 
satisfies: 

2RR 2 + 2R(l-f 2 )=F(r) } (3) 

with F being an arbitrary function of class C 2 , R = dR/dt and R' = dR/dr. We have three 
distinct solutions depending on whether f 2 < 1, = 1, > 1 and they correspond to elliptic 
(closed), parabolic (flat) and hyperbolic (open) cases, respectively. 

The proper density can be expressed as: 

F' 
P = 16tcR'R 2 ' ^ 

Whatever the curvature, the total mass within comoving radius r is: 

M(r) = - f r drj~ x F' = 4vr f drpf^R'R 2 , (5) 

4 jo Jo 

so that 

M '(r) = ^- = 47cpf- l R'R 2 . 
dr 

Also for p > everywhere we have F' > and R' > so that in the non-singular part of 

the model R > except for r = and F(r) is non-negative and monotonically increasing 
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for r > 0. This could be used to define the new radial coordinate r 3 = M(r) and find the 
parametric solutions for the rate of expansion. 

In the flat (parabolic) case f 2 = 1, we have 

R=^9F) 1/3 (t + pf 3 , (6) 

with (3(r) being an arbitrary function of class C 2 for all r. After the change of coordinates 
R(t, f) = f(t + /3(f)) , the metric becomes: 

ds * = df _ {t + pf/3 (y2 dr 2 + r 2^2^ ? (?) 

where 

Y = ! ■ 2 ^' 



3 (t + /?) ' w 

and from @ the density is given by 

P = ^—5-- (9) 

Clearly, we have that (t — > oo) the model tends to the flat Einstein-de Sitter case. 

For the closed and open cases the parametric solutions for the rate of expansion can be 
written as (Bonnor 1972; 1974): 

R = \f(\- f) l [1 - cosOO] , f < 1, (10a) 

t + (3 = ±F(l-f 2 y 3/2 iv-sm(v)], f 2 <l, (10b) 

^F(/ 2 -l) _1 [cosh(i;)-l], f 2 >l, (11a) 

l„/, 2 ^-3/2 , 2 



and 

t + /3=|F(/ 2 -lj "'"[sinh(i;)- V ], f > 1, (lib) 

with /5(r) being again a function of integration of class C 2 and v the parameter. 

The flat case (f 2 = 1) has been extensively studied elsewhere (Moffat and Tatarski 
1992). The model depends on one arbitrary function /3(r) and could be specified by 
assuming the density on some space-like hypersurface, say t = t . However, specifying the 
density on the past light cone of the observer is more appropriate. 

The cases of interest to us, ( |T0"D and (|TT|), correspond to closed and open models, 
respectively. 
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Before we proceed (in the next section) to discuss the observational grounds for 
modelling a local void, we need to amplify the discussion of the LTB model by introducing 
basic features of the propagation of light. The high degree of isotropy of the microwave 
background forces us to the conclusion that we must be located close to the spatial centre 
of the local LTB void. In our discussion, for the sake of simplicity, we place an observer at 
the centre (tob — to, r ob = 0). 

The luminosity distance between an observer at the origin of our coordinate system 
(£o,0) and the source at (t e ,r e , 9 e , </> e ) is (Bondi 1947): 

d L = (^=) = R(t e ,r e )[l + z(t e ,r e )} 2 , (12) 

where £ is the absolute luminosity of the source (the energy emitted per unit time in the 
source's rest frame), T is the measured flux (the energy per unit time per unit area as 
measured by the observer) and z(t e ,r e ) is the redshift (blueshift) for a light ray emitted at 
(t e ,r e ) and observed at (to,0). 

The light ray travelling inwards to the centre satisfies: 

ds 2 = dt 2 - R' 2 (t, r)f~ 2 dr 2 = 0, d6 = dtf> = 0, 

and thus 

^ = -R'(t,r)/f(r), (13) 

where the sign is determined by the fact that the light ray travels along the past light cone 
(i.e. if r e i > r e ", then t e > < t e "). 

Without getting into a detailed discussion, which can be found in Bondi 1947, or 
Moffat and Tatarski 1992, let us state that if the equation of the light ray travelling along 
the light cone is: 

t = T(r), (14) 

using (|i"3D , we get the equation of a ray along the path: 

dT(r) R' 



dr f 

where 



R'[T(r),r] °~ H 



dtdr 



[T(r),r], (15) 



dR' 



r,T(r) °t 



r,T(r) 

The equation for the redshift considered as a function of r along the light cone is: 

dz 

- = (l + ^'[T(r),r], (16) 
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and the shift z\ for a light ray travelling from (tx,ri) to (to,0) is: 

ln(l + zi ) = -ln(l-Oi - / dr — K J (17) 

jo r(l — aij 

where 

ai(r) = .R[T(r),r], 

and, in obtaining equation (|TTD, we used ([|) and (|5|). Thus we have two contributions 
to the redshift. The cosmological redshift due to expansion, described by the first term 
with a\ = R, and the gravitational shift due to the difference between the potential energy 
per unit mass at the source and at the observer. Obviously, in the homogeneous case 
(M'(r) = 0) there is no gravitational shift. 

3. LOCAL VOID 



If we restrict ourselves to spatial scales that have been well probed observationally, i.e. 
up to a few hundred Mpc, the most striking feature of the luminous matter distribution 
is the existence of large voids surrounded by sheet-like structures containing galaxies (e.g. 
Geller and Huchra 1989). These surveys (see also Efstathiou et al. 1990; Saunders et 
al. 1991) give a typical size of the voids of the order 50-60 h~ x Mpc. There has also 
been some evidence (Broadhurst et al. 1990) -with less certainty- for the existence of 
larger underdense regions with characteristic sizes of about 130 h^ 1 Mpc. Also, dynamical 
estimates of the FRW density parameter Qq give very different results on different scales. 
The observations of galactic halos on scales less than about 10 to 30 Mpc typically give (see 
e.g. Sancisi and van Albada 1987) fiio-30 — 0.2 ± 0.1. On the other hand, smoothing the 
observations over larger scales (> 20 Mpc, say ~ 100 Mpc ) indicates (e.g. Efstathiou et al. 
1990; Saunders et al. 1991) the existence of a less clustered component with a contribution 
exceeding 0.2, and perhaps as high as fi^ioo — 0.8 ± 0.2. 

At the same time, the large scale galaxy surveys (some of the recent literature is given 
in Maddox et al. 1990; Tyson 1988; Heydon-Dumbleton et al. 1989; Lilly 1993) firmly 
indicate a considerable excess in the number-magnitude counts for faint galaxies relative to 
predictions of homogeneous, "no-evolution" models. This excess could be the result of a 
non-standard galactic evolution or could be caused by rather exotic FRW cosmology (i.e. 
the deceleration parameter q <C 0.5 or a non-zero cosmological constant A). However, 
it can also be treated as an observational indication of a very large (on the scale of the 
redshift z ~ 0.5) void. In the following, we choose this latter option in interpreting the faint 
galaxies number counts and model the density distribution of a local void accordingly. (A 
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somewhat similar model of a gaussian void as a function of the comoving coordinate can be 
found in Moffat and Tatarski 1994.) 

We study a void with the central density equal to that of an FRW model with the 
density parameter Qq = 0.2, asymptotically approaching the FRW model with Qq — 1. Since 
cosmological observations are done by detecting some form of electromagnetic radiation, 
the mapping obtained from them describes the density along the light cone. We describe 
the density distribution, as a function of the redshift z, by: 



QJz) 



^rnin + {z/a) A 



a 



[z/aY 



(1* 



The choice of a rational function Q v (z) (as opposed to, say, a gaussian distribution) assures 
that the differential equations that we solve numerically (already quite complex) are not 
unnecessarily complicated further. In the numerical calculations presented in the next 



section we used the values fL 



0.2, a 



1. The two density distributions presented 



there are parametrized by (A) a = 0.125 and (B) a = 0.25, and are depicted in Figure [IJ. 
ll 




Fig. 1. — The density distributions Q v as functions of the redshift Z. The parametrization of the cases 
(A) and (B) is given in the text. 

In the following figures the results obtained for the (steeper) distribution (A) will be 
represented by dashed lines, while the results for the density distribution (B) will always be 
given as solid lines. Also, since we study here the observations of the luminous matter, we 



-7- 



present our results in the range of up to z ~ 6. (An observation of an object with a redshift 
in excess of this value would be hard to reconcile with standard FRW cosmology.) 

The final question we wish to address in this section is that of shell crossing in the LTB 
models. The problem has been extensively studied (for a detailed review and a complete 
list of references see Krasihski 1994). 

As seen from (|2|) and (|]) shell crossing can occur when R'(t,r) = for some r = r s . 
This, in general, may lead to lim r ^ rs p = oo and a shell crossing singularity. Different shells 
r = const collide and the comoving coordinates become inadmissible. The singularity can 
be avoided, if the functions F'(r) and f(r) both have at r = r s zeros of the same order as 
R'(t,r) has. This, however, is not the case in the present work (we assume F(r) > and 
F'{r) > 0). 

In the simple model presented here, shell crossing does occur. The physical significance 
of this can be viewed from diametrically opposite positions. One extreme is to regard a 
shell crossing singularity as unphysical and try to exclude it by assuming special initial 
conditions, as was the strategy in the work of Bonnor and Chamorro (Bonnor and Chamorro 
1990; 1991; Chamorro 1991). However, this leads to difficulties with adjusting the model to 
observational data. In short: either the average density of the void is not appreciably lower 
than the density of the background, or the matter in the void is much older than outside. 
The second approach, best represented by Sato and collaborators (see e.g. Sato 1984; Suto 
et al. 1984a, 1984b), is to treat the shell crossing singularity as a physical entity. In this 
approach, shell crossings -numerically modelled as surface layers of matter- are interpreted 
as rims of the expanding voids and correspond to the aforementioned sheet-like structures 
surrounding observed voids. 

An observationally acceptable LTB void, with the density distribution similar to that 
depicted in Figure [l], is too deep to escape shell crossing and, at the same time, allow for 
the simultaneous Big Bang throughout the whole space. Since we see the latter property 
as very important in order to be able to incorporate the accomplishments of the FRW 
cosmology in dealing with epochs preceding the matter dominated era, we choose to accept 
shell crossing as a price. The more so, that -as will be shown in the next section- the 
distribution (|TB"D, even though chosen on purely observational grounds, leaves us with a 
bonus: the shell crossing happens in the distant future. This, in our opinion, indicates that 
the LTB void studied here is applicable to the matter dominated era. (The fact that the 
LTB model does not allow for pressure is not perilous from our point of view. After all, the 
universe has been matter dominated since z ~ 10 4 . The model is applicable to this era and 
may be thought of as a continuation of an earlier -very nearly- FRW model, provided its 
density contrast at the time of decoupling, t = t dc , is p(t dc , r) / p FRW (td c ) ~ 1.) 



4. RESULTS 

In general, an LTB model depends on three arbitrary functions (see section |2Tf ), 
F(r),/3(r) and f(r). Since F(r) can be interpreted as twice the effective gravitational mass 
within comoving radius r (Bondi 1947), then, in accordance with the discussion following 
(|), assuming its form is equivalent to a coordinate choice. In our calculations we used 
F(r) = 4r 3 . The second function, /3(r), sets the initial singularity hypersurface of the 
model. Since we want the region of highest density in our model to be fully equivalent to 
the FRW universe, we set j3(r) = Po = const, thereby assuming a universally simultaneous 
Big Bang. In doing so, we give up a very important feature of an LTB model: an extra 
(with respect to FRW) degree of freedom that would allow the age of the universe to be 
position dependent. However, we study the following two cases (both with the simultaneous 
Big Bang hypersurface): the universe whose age is equal to that of the FRW critical 
(Qq = 1) case and the universe with the age equal to that of the FRW f2 = 0.2 one. The 
parametrization of the cases is, respectively, (I) (3q = Pi = and (II) (3q = j3n = 1.23606, 
where -in both cases- we set the time coordinate of constant time hypersurface "now" so 
that it is equal to the age of the universe to = 1 in the FRW model with Q = 1. The 
third ("curvature") function, f(r), is an unknown to be solved for in our calculations. Since 
we are modelling an underdense comoving void in an FRW universe, we choose the LTB 
hyperbolic (f 2 > 1) case. 

In a manner similar to that employed in Moffat and Tatarski 1992, we assume that since 
all cosmological observations are necessarily done by detecting some form of electromagnetic 
radiation, the solution should progress along the light cone. The final set of equations we 
solve consists of the equations QToD, flTED and the equation describing the ratio of the local 
density (Q) to the density of the FRW universe through the relation ( |l8l) taken along the 
past light cone: 

mw = 3 r'T5 r) +f o i 2 . = n. W r)]. (i9) 

PFBw[T(r)} SR'[T(r),r]R 2 [T(r),r] " [ K " ' 

The final form of the equations is quite unwieldy and we do not present them here. Since 
T{r) (the time of emission t e of a light ray observed at r = at to) is now given by (|llb|) 
and z(r) is governed by (|16D, the functions to be solved for are f(r),z(r) and v(r), where 
the parameter v becomes the function of position. 

The initial conditions for the integration have to be set at r ^ 0, since the analytic 
expressions (|ll|) are singular at r = 0, where f 2 = l (we have a flat Qo = 0.2 FRW universe 
there). We assume that for the initial radius r^ <^ 1 (we use dimensionless radius and time 
in the calculations) corresponding Zi and ti are given by their standard FRW values. Then 
z(ri),v(ri) and /(n) can be obtained from ( |llb| ), ( |T5| ) and (p!6|). 
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Once the equations have been numerically integrated (in all cases the accuracy of our 
numerical procedure was 10~ 6 or better □) we use ( |1 1 b| ) to obtain t e for a given r e . The 
luminosity distance di corresponding to this event is obtained with the use of (|T^) and 
z(t) (useful in studying the cosmological time scale) is given by the parametric relation 
[T(r),z(r)\. 

The results of our numerical calculations for z(r), where r is the dimensionless 
comoving radius used in the calculations are depicted in Figure 0. The coordinate distance 
has no direct physical relevance, but our units here are such that r = 1 corresponds to 
2997.95/i~ 1 Mpc, where h is the usual coefficient in the observationally determined value of 
the Hubble constant: H = lOOh km s _1 Mpc _1 . 




Fig. 2. — The redshift Z as a function of the comoving radius r. Roman numerals I and II denote the 
"younger" (/?o = /3j) and "older" (J3q = (3jj) solutions, respectively. Dashed and solid lines correspondence 
to the two density distributions of the void is the same as in Figure 1. FRW f2 = 0.2 and f2 = 1 solutions, 
respectively, are denoted "FRW 0.2" and "FRW 1". 

The behaviour of z(r) is as intuitively expected. In the "older" case (II) the redshift 
z asymptotically converges to the FRW Qq — 0.2 behaviour. The additional gravitational 



2 We do not elaborate here on the intricacies of the numerical procedure. A good discussion on the 
technical aspects of numerically solving a somewhat similar problem involving an LTB model can be found 
in Paczyhski and Piran 1990. 
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shift caused by the mass distribution of the void is apparent on the scales comparable 
to that of the void (z £ 0.5). On intermediate to large scales the effect of a cosmological 
expansion is predominant. In the "younger" case (I), z(r) -after an excess over the FRW 
relation due to the gravitational shift induced by the void- asymptotically tends to a limit 
that could, in accordance with the FRW interpretation, correspond to the universe with 
the density parameter in the range Q e (0.2; 1). In both instances, the convergence to the 
FRW relation is marginally faster for the cases with steeper density distributions. 

However, one should not forget that the comoving distance is not an observable, whereas 
the luminosity distance oIl is. In principle, provided we know its absolute luminosity £, we 
can establish the luminosity distance, defined by Andi 2 = C/J- ', by measuring the energy 
flux T of an observed object (for a discussion of usual caveats associated with so-called 
"standard candles" see e.g. Weinberg 1972; Kolb and Turner 1990; Peebles 1993). It is the 
z{d£) relation that is observationally relevant. Let us recall that the FRW equivalent of 

(0) is: 

d L 2 = a 2 (t )r e 2 (l + z e ) 2 . (20) 

Before we probe the observable properties of the solution, let us investigate the 
temporal evolution of the void. In doing so we seek answers to two important questions. 
First, how soon in the further evolution of the void the formation of a shell crossing 
singularity occurs. If this happens in the distant enough future, then our model is applicable 
to a sufficiently long cosmological epoch to be considered physical. Secondly, we want 
to verify that at early times the density contrast p v /pb ~ 1, where p v and pb are matter 
densities in the void and in the background, in accordance with our conjecture that the 
model be very similar to an FRW one at the beginning of the expansion. 

In general, the earlier the formation time of the void, the faster the growth of the 
density contrast. Also, the rate of growth is increased in the cases of higher initial contrasts. 
(See Sato 1984. A very thorough numerical study of the evolution of general relativistic 
voids, including the LTB ones, can also be found in Vadas 1993.) In the study presented 
here, however, we did not assume a density distribution on some hypersurface of constant 
time as initial condition. The observationally based density distributions of Figure [1] were 
given as functions of the redshift z and since we numerically integrated the equations along 
the light cone, we are not yet in the position to examine the time evolution of the spatial 
distribution of the matter density. However, equipped with the solution for /(r) we can use 
( |llb|) to obtain the values of the parameter v^ir) corresponding to a given hypersurface 



of constant time t = th throughout the space. Then, using (|J) we get the spatial matter 
density distributions for the time th- 

The results (numerical solution points) of this calculation for the future and past 
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evolution of the steeper void (A) in the "younger" case (I) are presented in Figures [3] and f|, 
respectively. The density distribution in this case evolves fastest and, as such, is the most 
representative in estimating the time scales involved. 

1, 




le-10 



100. R 1000. 



le+5 



Fig. 3. — The future evolution of the void (A) in the case (I): to + Pi = *0 = ^FRWX- The 
matter density as a function of R{thi r )/R(to, 7"j), where Tj is the initial radius, is shown at times 
t h = t , 10£ , 100£o and lOOOto- 

Figure |3] shows that the void expands outward increasing its density contrast and 
forming a very steep ridge on its boundary. This (in accordance with Sato 1984; Vadas 
1993) leads to the eventual formation of a thin dense shell and the shell crossing singularity. 
Of interest to us, however, is the time scale of this process. Since the shell crossing occurs 
for t > 1000io, where to is the present age of the universe, we conclude that our model is 
applicable to a satisfactorily long cosmological epoch. 

At the same time, Figure |] clearly indicates that our model satisfies the requirement of 
similarity to the FRW one at early cosmological times. At th = O.OOlto the matter density is 
approximately constant everywhere. For the case depicted in Figure f| the density contrast 
is p v /p b « 0.982 for t h = 0.001t (« 0.996 for t h = 10"%). In all cases at t/{t + p ) ~ 10~ 6 
we have \p v /pb — 1| — 10~ 5 . There is hardly a trace of the void's presence. 

The z(t) relations for the cases of different "ages of the universe" are presented 
in Figures ^| and |6|. There is some divergence from the FRW behaviour, but for early 
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1000. 

K 

Fig. 4. — The past evolution of the void (A) in the case (I): to + Pi = *0 = ^FRWl- The 
matter density as a function of R(th,r)/ i?(*o,rj), where 7"j is the initial radius, is shown at times 
t h = t , O.lto, O.Olto and 0.001£o- 



cosmological times the relations tend to their FRW counterparts. This asymptotic behaviour 
is in accordance with our assumption of a simultaneous big bang, (3(r) = const , and with 
our setting the age of the universe in each case to be equal to that of the corresponding 
FRW case. In both (I) and (II) cases objects with similar redshifts are younger than their 
FRW counterparts (i.e. t LTB > t FRW , when z LTB = zfrw), but not significantly. Those 
objects have marginally more time to form and evolve than their FRW counterparts. 

Let us now explore the z(di) relations resulting from our model. Any nonstandard 
cosmological model to be of interest must pass two tests when confronted with observations. 
To be considered realistic, it is of utmost importance that its predictions do not contradict 
the linearity of Hubble's law z = Hd, well established on small scales. At the same time, 
it should be sufficiently different from the standard (FRW) models to allow for significant 
reinterpretation of observational results. Let us now apply these two tests to our model. 

Figure [7] shows that (in all cases studied here) on large scales the LTB "Hubble 
relation" between the redshift z and the luminosity distance di is reminiscent of the FRW 
relation for the universe with the density parameter in the range Q G (0.2; 1). 

At the same time, inspection of Figure || shows that on small scales very nearly linear 
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'Hubble diagrams" are obtained for both types of LTB voids and both time scales. 
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z3 




Fig. 5. — The redshift Z as a function of the cosmic time t for the case (I): to + /3j 
corresponds to the to in the FRW critical case.) 



tFRWl- (t — 1 



z3 




Fig. 6. — The redshift Z as a function of the cosmic time t for the case (II): to + /3jj = tFRW0.2- (t — 1 
corresponds to the to in the FRW critical case.) 
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Fig. 7. — The Z vs. di relations for both the "younger" (I) and "older" (II) cases. The luminosity distance 
di, is given in Mpc. 
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00 200 400 600 d 800 1000 1200 1400 

Fig. 8. — The Z vs. di relations for both the "younger" (I) and "older" (II) cases for smaller scales. The 
luminosity distance d^ is given in Mpc. 
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If we recall (Sandage 1972, a good overview can be found in §5 of Peebles 1993) that 
the Hubble diagram, extended to z ~ 0.5, follows the Hubble's law with the redshifts' 
scatter (assuming identical absolute magnitudes of observed galaxies) of ~ 20 % around the 
theoretical prediction, then we can claim that the redshift-distance (redshift-magnitude) 
relations inferred from our model are within the observational bounds. Also, inspection of 
Figure |8| shows that on small scales observationally indistinguishable from linear "Hubble 
diagrams" are obtained. However, a different value for the Hubble parameter (constant) 
is inferred, if we insist on interpreting the results of cosmological observations through an 
FRW model. 

Also for larger observed redshifts (see e.g. the results of McCarthy's survey of the state 
of the measurements of the infrared Hubble diagram for radio galaxies, out to z ~ 4, in 
§5 of Peebles 1993) our results are well within the observational scatter. Moreover, when 
the measurements of the Hubble diagram are extended to large scales an interesting feature 
emerges. The slope of the diagram at z < 1 seems to be consistent with the linear Hubble 
law. However, at z > 1 a hint of curvature appears. In the standard FRW interpretation, 
this is presumed to be the result of relativistic corrections to the redshift-distance relation 
and of the fact that the galaxies observed at high redshifts are seen as they were younger 
and very likely more luminous than the galaxies observed at low redshifts. 

The fact that the redshift is observed to increase with distance, at least roughly in 
accordance with Hubble's law, for redshifts out to in excess of unity, is usually heralded 
as one of the major observational confirmations of standard FRW cosmology. It would be 
interesting to see how the predictions of our model withstand the same observational test. 

Figure [5] depicts the logarithmic redshift-luminosity distance relations for both LTB 
voids and both time scales. A survey of measurements (adapted from Lilly 1994) is also 
included. The z{d£) relations for the FRW cases with Oo = 1 and Qq — 0.2 are presented 
for comparison. 

It is now apparent that the redshift-luminosity distance relations inferred from our 
model are in agreement with observations. For small redshifts, say up to z ~ 0.1, they 
are observationally indistinguishable from the FRW results. For intermediate redshifts 
-in the range, say, z G (0.1; 1)-, where the gravitational shift contribution to the overall 
redshift causes the most manifest departure from the FRW relations, our results remain 
well within the observational scatter. In the range z > 1, where the "hint of curvature" in 
the redshift-distance relation is observed, our results fit the observations better than the 
FRW Q = 1 prediction. However, a still better fit is provided by the FRW Q = 0.2 curve. 

Due to the possible (e.g. evolutionary) corrections to the redshift-distance relation 
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Fig. 9. — The log(z) vs. \o<g[dj}j relations. Observations, denoted by O, are adapted from Lilly 1994. 
The luminosity distance d^ is given in Mpc. 



and the minute dissimilarities between FRW and LTB relations no assertions as to the 
best fit to observations seem irrefutable. Nevertheless, it is worth noticing that the model 
presented here supplies us with means to reinterpret the correlation between the theory and 
observational results (perhaps improving agreement between the two). 

As a final element of the discussion of observational properties of our model with a 
local LTB void, let us concentrate on the Hubble parameter (constant) measurement. In 
a spherically symmetric model we have two "Hubble parameters": H r (t,r) for the local 
expansion rate in the radial direction and H±(t,r) for expansion in the perpendicular 
direction. Usual definitions give: 

lr R , 

(21a) 



H r = — = — , 



H, 



(21b) 



l r R 

i± r 

h~R' 

where / denotes the proper distance, i.e. l r = R'(t,r)f~ 1 dr and l± = R(t,r)dQ. Due to the 
fact that there are both gravitational and expansion redshifts contributing to the total z, 
neither of the Hubble parameters H r , H^ is fully analogous to the FRW's H FRW = a/a. 
The closest analogy exists for small separations of the source and the observer. Let us 
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notice that the integral in ( W7h can be rewritten as: 



r i , A-npr 
dr 



o (l- ai ) 2 ' 

and for small T\ can be neglected when compared to the first term of fll7l) . Expanding the 
logarithms on both sides we see that light emitted at (t e ,r e ) and observed at (to,0) satisfies 
for small r e or small (to — t e ): 

z{t e ,r e ) = R(t e ,r e ), 

where t e is T(r e ) from (|15|) with the initial condition T(0) = to- Using ( |21b|) and fll2 ) we 
get for small r e : 

z(t e , r e ) = H ± (t e , r e )d L {t e , r e ), (22) 

which is formally analogous to the FRW result. Two main differences are that our relation 
is local and that from cosmological observations (on small scales) we obtain the angular 
Hubble parameter H± = R/R rather than Hfrw = a/a. 

In general, if we lived in a local LTB void and the z vs. d^ relation differed from the 
FRW one as described in this paper, but we were biased by our theoretical prejudice and 
interpreted cosmological observations through an FRW model, we would expect the value 
of the Hubble parameter to be position, or rather di,, dependent. 

To explore this possibility let us recall that in FRW cosmology the exact result for the 
Hubble relation (z versus di) in the matter dominated universe is (Weinberg 1972; Kolb 
and Turner 1990): 



Hodi = qo 



-2 



zq + (g - 1) ( \f2zq + 1 - l) 



(23) 



where q = —a(t )/a(to)H 2 is the deceleration parameter. 



On cosmologically very small distances we measure the same value of H independently 
of the model (we call this value "the local measurement"). This stems from the fact that, 
due to our assumptions, very close to the centre (r<l) the model is well approximated 
by the FRW universe with Q = 0.2. Obviously, if the universe were locally LTB rather 
than FRW, then the Hubble parameter based on the observed (LTB) values of z and d^ 
but inferred through an FRW relation (p3|), would be position (redshift) dependent. The 
dependence of the Hubble parameter H (in units of the H value as measured locally) on 



the redshift z is shown in Figure [10]. For the sake of the clarity of presentation, only the 
results for the "younger" case (I) are included in the figure. Qualitatively, the "older" case 
(II) exhibits the same behaviour. Quantitatively, as seen from the z vs. di relations of 
Figure |7| and (p3|), the "observed" values of the Hubble parameter (in the units of the local 
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Fig. 10. — The "observed" Hubble constant H (in units of the local measurement) as a function of the 
redshift Z. 



measurement) for the case (II) become increasingly smaller than those of the case (I) (~ 5% 
smaller for z«4), 

The H(z) relation depicted in Figure [H| strongly depends on the choice of the FRW 
case -the choice of the deceleration parameter q Q in ( p3|) ~ through which we interpret the 
observations. The main feature in the FRW £1 = 1 interpretation is a distinct monotonic 
correlation between H and z (or di,). The values of the "observed" Hubble constant decrease 
with the redshift z asymptotically approaching the background limit. If we interpret the 
results within the FRW Q < 1 framework, the "observed" values of the Hubble constant 
first decrease with z and then asymptotically increase to some background limit. (This 
asymptotic behaviour becomes clearly visible only for very large z.) The position of the 
minimum in H depends on the size of the LTB void. Neither of these patterns is evident in 
the reported measurements of the Hubble constant. 

However, one has to remember that the measurements of the Hubble constant come 
inseparable from their own observational uncertainties (estimates of the distance using a 
distance ladder, corrections for peculiar motion etc.). Moreover, it is not the broader idea 
of a local LTB void embedded in a globally FRW universe that fails here. A choice of 
density distributions Q(z) more complex than (18) could result in a better agreement with 
observations. 
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We think that the variance in H(z) exhibited in Figure 10 is a valuable feature of 
the model. The ratios of the highest to the lowest value of the "observed" H within the 
range of redshifts of Figure 110 are ~ 1.5. At the same time, the values of H Q reported to 



date span the range 40 to 100 km s _1 Mpc _1 (with standard errors quoted frequently as 
10 km s _1 Mpc _1 or less!). Inhomogeneities similar to the LTB void presented here might 
provide an explanation for this. 

5. DISCUSSION 



The LTB voids presented here decrease their density contrasts (the depth of the void 
with respect to the FRW background) when evolved back in time. At early times they are 
almost homogenized: at t/(t + /3q) — 10~ 6 we have \pLTB(r)/pFRw — 1| — 10~ 5 everywhere. 
This corresponds to a universe which at the beginning is very similar to the FRW one, but 
different at late stages. 

In a model utilizing a local LTB void, while retaining the accomplishments of the 
FRW cosmology in dealing with epochs preceding the matter dominated era, we can gain 
new freedom in modelling the more recent universe. We can solve the age of the universe 
problem by assuming (3{r) = const ^ 0, provide the excess power observed on scales of 
5 — 10,000 km s _1 in modelling structure formation (see Moffat and Tatarski 1992) and 
provide an explanation for the wide range of reported values of the Hubble constant. 

The use of a locally inhomogeneous model enriches the spectrum of available 
interpretations. The ultimate question to be addressed, however, is the agreement between 
predictions of the model and observations. As has been demonstrated by the results of 
the present work, the observationally based density distributions spanning the range of 
Q G (0.2; 1) satisfy the observational tests. The agreement with observations is achieved 
both in the case of the model with the "age of the universe" equal to that of the critical 
(Q = 1) FRW one and the "older" case with t + /3 equal to the FRW Q = 0.2 value. 
Theoretical prejudice (as powerful as the inflationary paradigm) might favour the critical 
FRW model. 
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was supported by the Natural Sciences and Engineering Research Council of Canada. One 
of the authors (DCT) thanks the Ministry of Colleges and Universities of the Province of 
Ontario for financial support through an Ontario Graduate Scholarship. 



-21- 



REFERENCES 

Bondi, H. 1947, M.N.R.A.S., 107, 410. 

Bonnor, W. B. 1972, M.N.R.A.S., 159, 261. 

Bonnor, W. B. 1974, M.N.R.A.S., 167, 55. 

Bonnor, W. B., and Chamorro, A. 1990, Ap. J., 361, 21. 

Bonnor, W. B., and Chamorro, A. 1991, Ap. J., 378, 461. 

Broadhurst, T. J., et al. 1990, Nature, 343, 726. 

Chamorro, A. 1991, Ap. J., 383, 51. 

Clowes, R. C, and Campusano, L. E. 1991, M.N.R.A.S., 249, 218. 

Efstathiou, C, et al. 1990, M.N.R.A.S., 247, 10. 

Geller, M. J., and Huchra, J. P. 1989, Science, 246, 897. 

Guth, A. H. 1981, Phys. Rev. D 23, 347. 

Heydon-Dumbleton, N. H., et al. 1989, M.N.R.A.S., 238, 379. 

Hogan, C. J. 1990, Nature, 344, 107. 

Kolb E. W. and Turner M. S. 1990, The Early Universe, Addison-Wesley Publ. Co., 
Redwood City. 

Krasihski, A. 1994, The Inhomogeneous Universe - Models and their Properties, to be 
published by Cambridge University Press. 

Lemaitre, G. 1931, M.N.R.A.S., 91, 490. 

Lemaitre, G. 1933, Ann. Soc. Sci. Bruxelles A53, 51. 

Lilly, S. J. 1993, Ap. J., 411, 501. 

Lilly, S. J. 1994, private communication. 

Lynden-Bell, D., et al. 1988, Ap. J., 326, 19. 

Maddox, S. J., et al. 1990, M.N.R.A.S., 240, 43. 

Mather, J. C, et al. 1990, Ap. J., 354, L37. 

Moffat, J. W., and Tatarski, D. C. 1992, Phys. Rev. D 45, 3512. 

Moffat, J. W., and Tatarski, D. C. 1994, to be published in the Proceedings of the XIV 
Moriond Workshop on Particle Astrophysics, Atomic Physics and Gravitation, 
Editions Frontieres, Paris. 



-22- 



Paczyhski, B., and Piran, T. 1990, Ap. J., 364, 341. 

Peebles, P. J. E. 1993, Principles of Physical Cosmology, Princeton University Press, 
Princeton. 

Sancisi, R., and van Albada, T. S. 1987, in Dark Matter in the Universe, eds. J. Kormendy 
and G. Knapp, Reidel, Dordrecht. 

Sandage, A. 1972, Ap. J., 173, 485. 

Sato, H. 1984, in General Relativity and Gravitation, eds. B. Bertotti, F. de Felice and A. 
Pascolini, Reidel, Dordrecht. 

Saunders, W. et al. 1991, Nature, 349, 32. 

Smoot, G. F. et al. 1993, Ap. J., 396, LI. 

Suto, Y., Sato, K. and Sato, H. 1984a, Progr. Theor. Phys. 71, 938. 

Suto, Y., Sato, K. and Sato, H. 1984b, Progr. Theor. Phys. 72, 1137. 

Tolman, R. C. 1934, Proc. Nat. Acad. Set., 20, 169. 

Tyson, J. A. 1988, A. J., 96, 1. 

Vadas, S. L. 1993, Phys. Rev. D 48, 4562. 

Weinberg S. 1972, Gravitation and Cosmology, John Wiley & Sons Inc., New York. 

Wright, E. L., et al. 1993, Ap. J., 396, L13. 



This preprint was prepared with the AAS IATj^X macros v3.0. 



